Plant Cell Reports

Improving regeneration of ribulose-1,5-bisphosphate (RUBP) is a promising approach to improve photosynthesis and plant growth. In addition to transgenic overexpression of target genes, it could be possible to directly overexpress endogenous target genes, through transcriptional enhancements. As shown by the recent discovery of a short sequence motif, that resembles the known octopine synthase (ocs) enhancer, transcriptional enhancement is achievable by relatively short endogenous sequences. In this study, we query the genome of several model and crop plant genomes for the presence of short enhancer motifs. We find hits across all genomes including some in promoter regions of genes. By using derivatives of these motifs in a transient fluorescence assay, we show that several of these are capable of inducing target gene expression in different promoter contexts. A motif scan of the created constructs, for the presence of known transcription factor binding sites, shows that the insertion of these motifs has created binding sites for different TGA-, NAC- and bZIP-transcription factors. Taken together our study shows the feasibility of finding enhancer sequences in the genomes of different plants. With advancement in gene-editing technologies, like prime editing, using such endogenous enhancer sequences, could allow for precise cisgenic promoter engineering of target genes.

Functional validation of genetic components in plants often requires cloning them separately into both plant and bacterial expression vectors, a process that is both time-consuming and laborious. This study aimed to simplify this workflow by developing plant-bacteria dual-host promoter systems that drive high-level constitutive expression in both environments. To achieve this, two variants of the chloramphenicol acetyltransferase promoter (PCAT), a bacterial {sigma} factor-dependent promoter, were integrated into the cauliflower mosaic virus 35S promoter (P35S), and their performance was evaluated using a hygromycin phosphotransferase (HPT)-GFP fusion reporter. One of these variants, PCAT1, conferred hygromycin resistance to Escherichia coli (DH5 and BL21 (DE3)) and maintained high-level expression comparable to the original P35S in onion epidermal cells. A hybrid P35S enhancer-PNOS system also conferred hygromycin resistance to E. coli, but its activity in inducing GFP signals in onion cells remained lower than that of P35S. Due to its compact size (89 bp) and efficiency, PCAT1 can serve as a module for converting standard plant vectors into dual-host systems, accelerating gene characterization and the development of new gene-based tools.

BackgroundMaize knobs are regions of constitutive heterochromatin that are readily identified in both meiotic and somatic chromosomes. These structures have been characterized as stable throughout the cell cycle, exhibiting late replication during the S-phase, and are composed of two specific families of highly repetitive DNA sequences: K180 and TR-1. Although widely used as cytogenetic markers due to their variability in number and chromosomal position across inbred lines, hybrids, and landraces, little is known about their chromatin structure and dynamics. In this study, we analyzed chromatin accessibility of knobs using DNS-seq data across four maize tissues representing distinct developmental stages. ResultsOur results reveal that K180 knobs exhibit tissue-specific variation in chromatin accessibility, transitioning between open and closed states during development. In contrast, the TR-1 knob of chromosome 4 remained consistently inaccessible across all tissues analyzed. A knob composed of both K180, and TR-1 further supported this observation, with only the K180 region showing dynamic accessibility. To validate these findings, we also analyzed other repetitive regions such as centromeres, which showed a uniformly closed chromatin structure similar to TR-1. These results suggest a unique developmental modulation of chromatin accessibility associated with K180 repeats. While the chromatin accessibility of knobs does not reach the levels observed at Transcription Start Sites (TSS), the comparison among different classes of repetitive DNA within maize constitutive heterochromatin provides compelling evidence for sequence-specific and tissue-specific chromatin dynamics. ConclusionsOur findings uncover a previously unrecognized property of maize knobs and establish a reference for future studies on chromatin organization and epigenetic regulation of repetitive DNA in plant genomes.

Triple-negative breast cancer (TNBC) remains one of the most aggressive breast cancer subtypes and lacks durable targeted therapies. Dysregulation of cell-cycle control, particularly through CDK4/6 signaling, is a defining feature of TNBC biology (Garrido-Castro et al., 2019). Extracts of Moringa oleifera have repeatedly been shown to induce G1-phase arrest in breast cancer models, yet the molecular basis of this phenotype remains unclear (Al-Asmari et al., 2015) (Gaffar et al., 2019). Emerging work on cross-kingdom regulation has raised the possibility that plant-derived microRNAs may, under specific conditions, interact with mammalian transcripts (Zhang et al., 2012) (Chin et al., 2016). Sequence shuffling for the negative control was performed with set.seed(42) to ensure reproducibility. Additional visualisations (nucleotide alignment and thermodynamic analyses) were generated using Python 3 (matplotlib v3.7). Here, we performed a high-stringency computational screen of conserved Moringa microRNAs against 30 genes implicated in TNBC pathogenesis using local sequence alignment. We identify a predicted high-affinity interaction between mol-miR156 and the human CDK4 3' untranslated region (3'UTR), characterized by an uninterrupted 12-nucleotide complementary motif that exceeds canonical mammalian microRNA seed requirements. These findings support the hypothesis that conserved plant microRNAs may exhibit latent structural compatibility with oncogenic human transcripts. While physiological delivery and functional repression are not demonstrated here, this work establishes a molecular framework for future experimental investigation into cross-kingdom RNA interactions relevant to cancer cell-cycle regulation. Impact StatementA high-stringency computational screen identifies latent molecular compatibility between a conserved plant microRNA and the human CDK4 oncogene, establishing a testable framework for cross-kingdom RNA interference in triple-negative breast cancer.

Plant hormones regulate almost every aspect of plant growth and development. KARRIKIN INSENSITIVE 2 (KAI2)-dependent signaling, which is thought to transduce signals derived from an unidentified ligand known as the KAI2 ligand (KL), regulates numerous traits, including seed germination in angiosperms and vegetative reproduction in bryophytes. The origin of KAI2 is believed to be ancient, and the evolution of its signaling pathways remains of significant interest. Ferns represent critical lineages for elucidating the evolution of land plant traits and growth mechanisms that enabled adaptation to terrestrial environments. Therefore, functional studies of key components of this pathway in ferns are essential for understanding the evolutionary trajectory of KAI2-dependent signaling during vascular plant diversification. However, experimental platforms for the CRISPR/Cas9 system, a powerful tool for investigating gene function, remain undeveloped in ferns. Here, we report an efficient CRISPR/Cas9 system based on ribozyme-gRNA-ribozyme (RGR) technology in the model fern, Ceratopteris richardii (C. richardii). We generated loss-of-function mutants of the KAI2 ortholog in C. richardii (CrKAI2), as well as the signaling components CrMAX2 and CrSMXL. We demonstrate that exogenous application of an artificial KL agonist increases the expression of KAI2-dependent signaling responsive genes in wild type plants; this response is abolished in Crkai2 mutants. These findings indicate that KAI2-dependent signaling is conserved in C. richardii. Furthermore, this study proposes an efficient CRISPR/Cas9 method that will facilitate genetic studies in ferns.

Long non-coding RNAs (lncRNAs) have emerged as critical players in plant development and stress responses, yet their involvement in pollination responses is largely unknown. To address this gap, we identified and characterized lncRNAs and their cis-acting, trans-acting, and miRNA-mediated regulatory interactions during both compatible and incompatible pollination in Arabidopsis thaliana. Leveraging publicly available datasets, we analyzed expression profiles at 10 and 60 minutes post-pollination. We identified 1,073 novel and 3,422 annotated lncRNAs, with 1,002 novel and 985 annotated, respectively, showing detectable expression after filtering. Differential expression analysis identified 12 lncRNAs at 10 min and 32 lncRNAs at 60 min post-pollination. Further investigation revealed 9 cis-targets, 112 trans-targets, and 144 miRNA-mediated regulatory interactions, many of which were enriched in pathways related to stress, defense, and self-incompatibility. Notably, the regulatory landscape is more active at 60 minutes than at 10 minutes post-pollination. These findings provide a robust framework and resource to facilitate future functional studies of lncRNAs during pollination.

O_LISymbiosis between legumes and rhizobia is beneficial on nutrient-poor soils, as it enables the fixation of atmospheric N2. To establish this symbiosis, gene expression in both the host plant and the symbiont has to be regulated. To understand the underlying RNA-mediated regulation of host gene expression, we designed experiments to identify competing endogenous networks involving circular RNA, microRNA, and linear transcripts during symbiosis, using wt and symbiosis-deficient Lotus japonicus mutants with the rhizobium Mesorhizobium loti (M. loti). C_LIO_LICircRNA, miRNA, and linear transcripts were identified from Lotus japonicus wildtype and CCamK mutant (ccamk-13; snf-1) seedlings without inoculation or with M. loti inoculation using deep short-read sequencing with rRNA-depletion and random primers. C_LIO_LIDifferentially expressed miRNAs showed negative correlations to predicted target genes and may regulate symbiotic processes. The symbiosis essential iron-sensor LjnsRING/BRUTUS expresses a circRNA which was upregulated in symbiotic treatments. This circRNA may act as a target mimic and contribute to nodule longevity. CircRNAs are predicted to act predominantly as trans-regulatory molecules with similar frequencies in Arabidopsis thaliania, Oryza sativa, and Lotus japonicus. C_LIO_LIWe identified novel miRNAs, long noncoding RNAs, and circRNAs, and nominated several as potential new regulatory non-coding RNAs that may act as target mimics to stabilize genes and support symbiosis. C_LI SummarySymbiosis between Lotus japonicus and Mesorhizobium loti involves treatment-specific regulation of competing endogenous RNA networks involving circular RNA, miRNA, and linear transcripts.

Understanding the regulation of enzyme activity involved in photosynthesis is essential for engineering enhanced carbon fixation in crops. In C4 plants, the enzyme phosphoenolpyruvate carboxylase (PEPC, EC 4.1.1.31) is one of the most abundant leaf enzymes and plays an essential role in photosynthetic carbon dioxide (CO2) fixation. The enzyme also plays a key role in central metabolism (e.g., providing intermediates to the citric acid cycle) and therefore must be highly regulated to coordinate its activity. The regulation of PEPC activity can occur allosterically by glucose 6-phosphate activation and malate inhibition, which is in part influenced by reversible phosphorylation. A specific light-dependent phosphorylation of PEPC at an N-terminal serine residue by the PEPC-protein kinase (PEPC-PK) can regulate its sensitivity to this allosteric regulation. However, the impact of this PEPC phosphorylation has not been tested in a C4 crop. Therefore, we created PEPC-PK mutant lines in Zea mays to assess the impact of PEPC phosphorylation on its allosteric regulation, photosynthesis, and growth. While the maximum PEPC activity was unchanged, PEPC in the PEPC-PK mutant plants was not phosphorylated under light and was more sensitive to malate inhibition. However, gas exchange, electron transport, and field biomass analyses showed no differences in the PEPC-PK mutant plants. These results demonstrate that in Z. mays PEPC phosphorylation affects enzyme sensitivity to malate in vitro but does not substantially alert photosynthetic performance or growth under field conditions suggesting additional regulation of PEPC activity in planta.

Hemp (Cannabis sativa) produces a wide array of medicinally significant compounds, including cannabidiol (CBD). These compounds are predominantly synthesized in female hemp inflorescences. The proposed research utilizes next-generation sequencing-based transcriptome analysis using a 3{square}-end-directed approach to identify differentially expressed genes between male and female hemp plants at the early vegetative stage. 886 differentially expressed genes (DEGs) were identified, a majority of which were upregulated in males compared to females. We hypothesized that alternative RNA processing contributes to sex-specific gene expression. To this end, 932 genes were identified that exhibited significant changes in poly(A) site usage when comparing males and females. These genes were much more likely to be differentially expressed, supportive of this hypothesis. Males tend to have longer 3 UTRs with canonical motifs found in the Near-Upstream Elements (NUE), compared to the shorter 3 UTRs in females, which have A-rich motifs near the cleavage site. This suggests that polyadenylation remodels hemp mRNAs with distal poly(A) sites being preferred in males. To further investigate when this sex-specific gene expression program is established, RNA was isolated from plants at various developmental stages, such as developing seeds, four-day-old seedlings, and different developmental stages up to four weeks after sowing. Diagnostic male-specific genes were analyzed using RT/PCR. The results indicate that sex-specific gene expression is not evident in seeds but rather is set during or after germination. SignificanceO_LIHemp males tend to have longer 3 UTRs with canonical motifs found in the Near-Upstream Elements (NUE), compared to the shorter 3 UTRs in females, which have A-rich motifs near the cleavage site. C_LIO_LIThe sex-specific gene expression program is not yet established in mature seed but is set in the time between germination and 4 days of growth. C_LI

Calreticulins are multifunctional proteins involved in calcium homeostasis, protein folding, and cellular signaling. In common bean (Phaseolus vulgaris), the molecular mechanisms that regulate infection and nodule development remain incompletely understood. The main objective of this study was to characterize the role of the calreticulin gene PvCRT08 during infection and nodulation processes. We first analyzed the calreticulin gene family in the P. vulgaris genome and identified three members, with PvCRT08 showing the highest transcript accumulation in roots and after inoculation with rhizobia. Spatial and temporal promoter analyses in transgenic composite bean roots revealed PvCRT08 activity in root hairs and in infected cells and vascular bundles of mature nodules. RNA interference (RNAi)-mediated PvCRT08 down-regulation in transgenic roots increased the number of infection threads and enhanced nitrogen fixation efficiency, leading to the formation of larger and more functional nodules, although total nodule number was unaffected. In contrast, overexpression of PvCRT08 impaired infection thread progression, reduced the expression of key nodulation marker genes (PvCyclin and PvNIN), decreased nodule number, and diminished nitrogen fixation capacity. These findings identify PvCRT08 as a key regulatory component of early infection events and nodule development in common bean. Furthermore, the study provides new insights into the molecular control of symbiotic efficiency and highlights PvCRT08 expression is critical to optimize the equilibrium between infection efficiency and nodule functionality.

RNA interference (RNAi) shows great potential to protect crops against fungal diseases, yet reported protection efficiencies vary greatly, and our understanding of the factors responsible for this variance remains limited. In this meta-analysis, we evaluated 89 studies that compare the efficiency of host-induced gene silencing (HIGS) and spray-induced gene silencing (SIGS) in controlling fungal diseases, focusing on biotrophic, hemibiotrophic, and necrotrophic fungi, the use of formulations, and the dsRNA design as explanatory factors for differences between reported efficiency values. Our results indicate that SIGS is slightly more effective, particularly in biotrophs. Surprisingly, SIGS studies using formulations did not outperform those applying naked dsRNA. We also assessed parameters of RNA design. Differences in dsRNA length and the number of constructs, and number of targets showed no consistent significant effect on resistance in either HIGS or SIGS. Interestingly, however, HIGS studies reported significantly higher efficiency when targeting genes closer to the 3 end and SIGS when targeting genes closer to the 5 end. We discuss potential reasons for the reported patterns, such as variability in dsRNA uptake mechanisms, intercellular trafficking and Dicer processing, and conclude that more research is needed to understand the biological mechanisms determining RNAi efficiency for fungal control.

The model plant species Arabidopsis thaliana has two nucleolus organizer regions (NORs), each consisting of tandem repeats of 45S rRNA genes on chromosomes 2 (NOR2) and 4 (NOR4). In the ecotype Col-0, the rRNA gene subtypes mapping to NOR2 are mostly silenced during development, whereas those mapping to NOR4 are mostly active. Using molecular and genetic mapping approaches, we demonstrate in multiple epigenetic mutants the occurrence of NOR conversions involving loss of rDNA and the associated telomeres in one NOR and replacement of the lost sequences with the corresponding sequences from the other NOR. We studied mutants of evolutionarily conserved chromatin remodelers DECREASE IN DNA METHYLATION 1 (DDM1) and CHROMATIN ASSEMBLY FACTOR 1 (CAF-1), and plant-specific CHROMOMETHYLASE 2 (CMT2) DNA methyltransferase. We observed NOR conversions in ddm1 and caf-1 mutants, where such NOR conversions altered the rRNA variant expression patterns, which often confound the data pertaining to release of rRNA gene silencing. We delineated the effects of these mutations on the rDNA instability-mediated alterations in rRNA variant expression from their effects on the actual release of rRNA gene silencing. We also show that these mutations release rRNA gene silencing independently of their effect on rDNA genomic instability involving NOR conversions.

Low temperature stress causes significant damage to the strawberry plant. During cold stress, plants undergo morphological and physiological changes often regulated at the genetic and/or epigenetic levels. Some strawberry cultivars are more cold-hardy than others. Using the diploid woodland strawberry as a model, we analyzed the effects of cold acclimation on methylome and transcriptome dynamics in the crowns and leaves of three ecotypes with contrasting cold tolerance. Alta, which was the most cold-tolerant ecotype, exhibited the highest genetic and epigenetic plasticity in response to cold. CHH-context methylation dominated the differentially methylated regions (DMRs) with more hypomethylation in crowns and hypermethylation in leaves. CG methylation was enriched in gene bodies, while non-CG methylation was prevalent in upstream and downstream regions. Our study revealed that less than a quarter of differentially methylated genes (DMGs) showed changes in transcript accumulation levels. This finding indicates that universal cold response in Fragaria vesca, as reflected by gene expression, cannot be mechanistically attributed to DNA methylation. The majority of differentially expressed differentially methylated genes (DEDMGs) were ecotype- and tissue-specific. Enrichment analysis revealed that these genes were involved in pathways related to stress tolerance, such as carbohydrate metabolism, lipid metabolism, ATP hydrolysis, and cellular detoxification. Each ecotype responded to cold through mobilization of its own set of differentially expressed genes (DEGs), DMGs, and DEDMGs, and variation in expression and methylation patterns exhibited by Alta, FDP817, and NCGR1363 suggest that cold signaling processes and survival depend on the tissue, ecotype, and geographical origin of the plants exposed to cold stress. Therefore, this study highlights the potential of both genetic markers and epialleles as molecular markers for the development of cold-tolerant octoploid strawberry cultivars that are better suited for propagation in Nordic climates.

BackgroundRAB Guanine Nucleotide Dissociation Inhibitors (RAB GDIs) are important vesicle transport regulators in eukaryotes, participating in the functional cycle of RAB GTPases by stabilizing their non-active GDP-conformation. AimsWe address the importance of the three Arabidopsis thaliana RAB GDI paralogs by genetic and developmental analyses and put these results into the seed plants evolution context. MethodsWe use methods of genetics, microscopy and phylogenetics. ResultsOur genetic analyses of Arabidopsis T-DNA insertional mutants confirm recent CRISPR alleles data indicating lethality of double gdi1 gdi2 mutants, and our microscopic data point to embryo development arrest in double mutant seeds. We also confirm the involvement of GDI2 and GDI3 in pollen tube growth. Moreover, our data show that GDI1 also contributes to proper pollen function. Our phylogenetic analysis reveals independent diversification of RAB GDIs in Gymnosperms and Angiosperms, with early specialization of an Angiosperm reproduction-and gametophyte-related clade. ConclusionsIn Arabidopsis, RAB GDI1 and 2 are important for the vegetative growth while RAB GDI2 and 3 are vital for reproduction. Evolution of the RAB GDI family reflects the evolution of seed plants. HighlightsRAB GDIs are vital for plant growth and reproduction and act redundantly. Even the low-transcribed RAB GDI1 isoform contributes to the proper pollen function. Two RAB GDI clades evolved in early Angiosperms.

Improving rice (Oryza sativa L.) yield requires a balanced enhancement of both sink size and source capacity. While many QTLs for sink size have been identified, only a few are known for source capacity, which is essential for achieving high yield. Here we identified qHP10 as a major QTL for increased photosynthetic rate by using chromosome segment substitution lines derived from a cross between the high-yielding indica cultivar Takanari and the average-yielding japonica cultivar Koshihikari. High-resolution mapping combined with CRISPR/Cas9-induced mutagenesis revealed that the causative gene underlying qHP10 is Mitogen-Activated Protein Kinase 4 (OsMPK4). A near-isogenic line carrying the OsMPK4Takanari allele (NIL-OsMPK4) had a 15-25% higher photosynthetic rate than Koshihikari. NIL-OsMPK4 also had higher stomatal conductance than Koshihikari but similar stomatal pore size and density, indicating that increased stomatal aperture increases photosynthetic rate. This enhancement is likely attributable to the down-regulation of OsMPK4 expression, which increases stomatal conductance and thus promotes CO2 uptake. Our findings demonstrate that OsMPK4 is a promising genetic target for increasing source capacity and, potentially, rice yield through molecular breeding. (175 words)

Particle bombardment systems are widely used for plant transformation, but commercial devices are expensive and rely on high-pressure helium gas. This study aimed to develop a cost-effective and helium gas-free alternative using an air duster gun connected to a commercial compressor. A nozzle (for DNA with transgenes), gold particles (as DNA carriers), nozzle-to-sample distance, and a method for coating gold particles with DNA were optimized to yield better transformation efficiency in targeting onion epidermal cells and rice calli. From the rice calli transformed with the newly developed system (a tool to shoot genes with massive air from a compressor: TSGMAC), stable transgenic plants could be obtained. TSGMAC offers a low-cost and helium gas-free solution for plant transformation and genome editing and can enhance accessibility to particle bombardment-based techniques.

Expression of the chloroplast AAA+ chaperone CLPD gene increases during senescence and drought, but its functional role in chloroplast proteostasis is poorly understood. This study provides a comprehensive analysis of Arabidopsis CLPD protein accumulation across development from early seedlings to senescence, and compares results to its homologs CLPC1,2, as well as CLPB3 and cpHSP90. The developmental consequences of complete loss of CLPD expression (clpd-1), as well as overexpression of functional CLPD or CLPD impaired in ATP hydrolysis (CLPD-TRAP), were determined in Arabidopsis. clpd-1 has accelerated seedling development while functional CLPD overexpression lines, but not CLPD-TRAP, have delayed development. To determine if CLPD is a bona fide CLP chaperone associating with the CLPPRT protease and to identify in vivo candidate substrates, we employed the CLPD-TRAP line during the vegetative and flowering (senescent) growth stages. Affinity purification of CLPD-TRAP followed by mass spectrometry showed high enrichment of the CLP protease complex, thus providing direct support for the role of CLPD in substrate delivery to the CLP protease. CLPC1,2 were also highly enriched in the CLPD-TRAP interactome, suggesting hetero-oligomerization and cooperation between the three chaperones is likely. Nine chloroplast candidate substrates were identified in the CLPD-interactomes, including: FHY2 involved in riboflavin synthesis, THI1 and THIC involved in thiamin metabolism, and four proteins of unknown function. Several of these have been previously identified as potential CLPC1 substrates; however, others appear to be specific to CLPD. CLPD acts in substrate selection within a heteromeric CLPC-CLPD hexamer, likely to make unique contributions through its divergent N-terminus.

Short Interrupted Repeats Cassettes (SIRC) are recently discovered eukaryotic DNA elements possessing many traits of satellite DNA and mobile genetic elements, and consisted of short direct repeats interspersed with diverse spacer sequences. The SIRC ensemble of individual species is highly heterogenous and cannot be studied using alignment methods. It was found that number of similar SIRC sequences in a given pair of species is in general correlated with their taxonomic distance, and, at the same time, closely related species can possess very diverged SIRC ensembles, which makes SIRC evolutionary pattern closer to mobile genetic element type. The SIRC sequences make up clusters with comparable sequence patterns, that are likely to demonstrate doublet evolutionary model which strongly supports that the SIRC structure is supported by the evolutionary selection. Several SIRC sequences of Arabidopsis were found to be of ancient origin with traceable evolution history as far as to the moss clade. We carried out unbiased detection of SIRC ensembles in 10 plant genomes and found that, despite very high intraspecies heterogeneity, SIRC sets possess strong interspecies phylogenetic signal. Key messageShort Interrupted Repeats Cassettes are elements of ancient origin, and could potentially be used to trace organism history, and to facilitate syntheny and Hi-C analysis.

Biolistic transformation is a versatile tool in plant science, yet high equipment costs and tissue damage from high-pressure gas remain significant barriers. Building on our previously developed "TSGMAC", a low-cost, helium-free biolistic system, we report three major advancements to enhance its throughput, delivery quality, and quantitative capability. First, a "guide barrel" assembled from commercial DIY fittings was developed; it effectively eliminates physical tissue damage and ensures uniform particle distribution, even in soft tissues like bok choy (Brassica rapa subsp. chinensis). Second, a rapid gene expression platform using PCR products was characterized. Results demonstrate that linear DNA constructs are efficiently circularized via non-homologous end joining (NHEJ) in plant cells, and protein expression is robust regardless of the relative positions of the promoter, coding sequence, and terminator. This system bypasses time-consuming cloning. Third, a cost-effective, highly sensitive dual-luciferase assay system utilizing teal Luc (teLuc) and inexpensive firefly luciferase (FLuc) inhibitors was established. This integrated workflow enables rapid, quantitative molecular biology using supermarket-obtained materials and standard PCR reagents. Our findings provide a practical foundation for plant scientists, synergistically accelerating gene functional analysis and genetic tool development.

(1) BackgroundSaffron crocus (Crocus sativus) is the source of saffron, the most expensive spice in the world. It evolved about 3000 years ago as a sterile triploid clone in Greece. Since then, saffron has spread across the globe, where regionally distinct practices of saffron cultivation have developed. Despite differences in morpho-physiological traits, genetic variability is low, if present at all. Here, we aim to resolve chromosomal and sequence-associated variability across saffron crocus cultivars from the crops main cultivation areas in Africa, Asia and Europe. (2) MethodsWe used genome-wide DNA polymorphisms obtained through genotyping-by-sequencing (GBS) of 33 saffron and 14 closely related Crocus accessions, which we place into a phylogenetic context. For karyotyping, we compare nine saffron accessions by multi-color fluorescent in situ hybridisation (FISH) with repetitive DNA probes. (3) Key resultsPhylogenetic analyses confirmed the single origin and clonal nature of all saffron accessions. We detected slight DNA differences among saffron crocus genotypes, which were minor compared with those in wild C. cartwrightianus populations. Still, the Iranian saffron accessions form a genetically very narrow group that differs from the other proveniences in population genetic analyses. However, chromosomes of some saffron accessions display variable FISH signals, likely resulting from gains and losses of tandemly repeated DNA. (4) Main conclusionsBased on the high genetic identity and small karyotypic differences, we confirm the clonal origin of the saffron accessions. Nevertheless, as we detected small and regional chromosomal variability, we conclude that at least four somaclonal saffron lineages emerged after saffrons origin. Societal Impact StatementFor millennia, many cultures developed cultivation practices and regional crop varieties. A notable case is saffron, the worlds most expensive spice that is harvested from stigmas of saffron crocus. This flower crop arose 3000 years ago in a singular genome triplication event and since then spread clonally across the globe. By identifying genetic and chromosomal variability in clonal saffron accessions, we highlight regional diversity, support the preservation of traditional knowledge, and underscore the risk of relying on only one clonal lineage. This informs strategies for saffron cultivation, linking cultural heritage with modern genomics to address biodiversity, evolution, and food security.